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The model of spatial-temporal distribution of point defects in a three-layer stressed nanoheterosystem 
GaAs/In x Gai_ x As/GaAs consideringthe self-assembled deformation-diffusion interaction is constructed. Within 
the framework of this model, the profile of spatial-temporal distribution of vacancies (interstitial atoms) in the 
stressed nanoheterosystem GaAs/In x Gai_ x As/GaAs is calculated. It is shown that in the case of a stationary 
state (f > 5t®), the concentration of vacancies in the inhomogeneously compressed interlayer is smaller rela¬ 
tive to the initial average value AT® by 16%. 
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1. Introduction 

Intensive development of nanotechnologies has provided an opportunity to create nanoelectronic 
devices on the basis of stressed nanoheterosystems GaAs/In x Gai_ x As/GaAs (ZnTe/Zn | _ x Cd x Te/ZnTe). The 
active region of such structures are layers In x Gai_ x As, Zni_ x Cd x Te, in which the electron-hole gas is 
localized being bounded on two sides of the potential barriers GaAs (ZnTe). 

It is known that optical and electric properties of such devices depend significantly on both the lattice 
deformation of the contacting systems and the spatial distribution of point defects. 

Such defects can penetrate from the surface or arise in the process of epitaxial growth. Besides, dif¬ 
fusion processes play an important role in the technology of fabricating optoelectronic devices. They are 
related with the redistribution of impurities in a semiconductor structure caused by both the ordinary 
gradient concentration of defects and the gradient of deformation tensor. 

The interaction of defects with the deformation field, created by both the mismatch of the crystal lat¬ 
tice of the contacting materials and the point defects, causes a spatial redistribution of the latter. It can 
lead both to accumulation and to a decrease of the number of defects in the active region (In x Gai_ x As, 
Zni_ x Cd x Te) of the operating element depending on the character of the deformation created both by 
the mismatch between parameters of contacting crystal lattices (co = 7% (4%) for GaAs/In x Gai_ x As/GaAs 
(ZnTe/Zni_ x Cd x Te/ZnTe), respectively (l]|2]), and by the action of defects. In particular, it is known that 
the gallium arsenide grown using the method of the molecular-beam epitaxy at low temperature contains 
an excess of arsenic 01 SI- Introduction of excessive arsenic causes a tetragonal distortion of the lattice 
material GaAs and the generation of point defects in it: interstitial atoms (As), vacancies (Ga) and anti- 
structural defects (Asoa), which, in turn, leads to their spatial redistribution. The lattice deformation and 
concentration of point defects generated under the action of the gradient of deformation tensor depend 
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on the mismatch between the lattice parameters of contacting layers of heterostructure, the temperature 
of growth, molecular fluxes Ga and As, concentration and chemical nature of the doped impurities. 

The strain caused by the mismatch between the lattice of the epitaxial layer and the substrate can 
be elastic when the thickness of the layer does not exceed a defined critical value 0- Otherwise, mis¬ 
match dislocations are formed accompanied by a sharp worsening of both the optical and the electric 
characteristics of devices. However, in the layers In x Gai_^As with the mismatch less than critical there 
is a significant decline of the mobility and the intensity of photoluminescence at certain terms 0, which 
is related to the increased number of point defects and a corresponding increase of the diffusion barrier 
to the atoms of the third group. 

In experimental work (6jJ, it is shown that in a heterostructure GaAs/In x Ga|_ r As, the stressed 
quantum-size heterolayers hamper the diffusion of hydrogen and defects into the bulk of the material 
which leads to a substantial difference of their spatial distribution in a heterostructure and homoge¬ 
neous layers. Theoretical research of the stationary distribution of defects within the framework of the 
self-assembled deformation-diffusion model has been considered in the work | 7]. 

Therefore, in order to create devices with prescribed physical properties, it is necessary to construct 
a spatial-temporal deformation-diffusion model that describes the self-assembled deformation-diffusion 
processes in stressed nanoheterostructures having their own point defects and impurities. 

The aim of this work is to construct a spatial-temporal deformation-diffusion model and calculate 
the spatial-temporal profile distribution of point defects (interstitial atoms and vacancies) in three-layer 
stressed nanoheterosystems GaAs/In, Ga | _j As/GaAs (ZnTe/Zn | _ A Cd,Te/ZnTe). 

2. The model of spatial-temporal redistribution of defects in a three- 
layer stressed nanoheterosystem 

Let us consider stressed nanoheterosystems GaAs/In,-Ga | _ ,.As/GaAs (ZnTe/Zn, _ x Cd Y Te/ZnTe) having 
interlayers InAs (CdTe) of the thickness 2a, that include three layers (figure 111, where N^, are 

the initial average defect concentrations, respectively, and D i, D 2 , />>, are diffusion coefficients. Suppose 
that the external layers GaAs (ZnTe) are of the thickness that considerably exceeds the width of the in¬ 
terlayer of the heterostructure (2 alL « 1), so deformation of these layers can be neglected [e,- (z) = 0, 
i = 1,3]. 

The mechanical deformation that occurs due to the mismatch between the lattice parameters of con¬ 
tacting materials of a heterosystem is approximated by the function (7): 



(1) 


where £0 = e xx + £yy + £zz < 0 is the relative change of the elementary cell volume of the grown layer on 
heteroboundaries z = |a|; e yy - e zz = (fl ,- + 1 - af)! a;, e xx = — (2C® /Cj^)£ yy , where i - 1,3 corresponds 
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Figure 1 . The model of stressed nanoheterostructures GaAs/In^-Gai-^As/GaAs (ZnTe/Zni_ x Cd 5: Te/ZnTe). 
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to the layers GaAs (ZnTe), i = 2 corresponds to the interlayers In,-Ga|_ x As (Zn|_ x Cd,Te); a,- are the crys¬ 
tal lattice parameters of the contacting materials GaAs (ZnTe) and InAs (CdTe) of the heterostructure, 
respectively; C® and C^' are the elastic constants of the material In x Gai_ x -As (Zni_ x Cd x Te). 

Epitaxial growth on a substrate with the mismatch between the lattice parameters takes place si¬ 
multaneously with the diffusion process, which is caused by both the concentration gradient of point 
defects [gradAT^ [z, t)] and the gradient parameter of deformation [gradf/,- (z, f)]. The latter induces an 
additional diffusion flux of defects, which is opposite to the ordinary gradient concentration flux of de¬ 
fects. Therefore, in the basis of this model it is necessary to put a self-assembled system of non-stationary 
equations for the parameter deformation Uj (z , t) and concentration of impurities [z, t) in a stressed 
heterosystem, the redistribution of which is performed similarly to an ordinary diffusion flux. 


t) = -Di 


dN™ (z, t) 
dz 


thus, by the deformation component of the flux 


C(z,t) = ~ Dl 


hr 


a dz 


where 6 l2) = K il) AQ W is the mechanical deformation potential, K M = (C{j + 20^) 13 is the module of 
uniform compression of the z'-th material, AQ W is the variation of elementary cell volume at the presence 
of a defect in the z'-th layer. 

Let the point defects be distributed with the initial average concentration nQ in the z'-th layer in a 
particular heterosystem. As a result of their self-assembled interaction through the deformation field, cre¬ 
ated by both the mismatch between lattice parameters of contacting materials of a heterosystem and the 
presence of defects, there is a variation of the concentration profile of point defects and of the character 
of deformation. 

The mechanical stress in epitaxial layers created by both the point defects and the mismatch between 
lattice parameters of contacting materials is described by the expression: 

Oi(z, t) = picjlJiiz, t) - 9^N [ J\z, t ) - pic^Eiiz), (2) 


where p;, a are the density of the z'-th medium and the longitudinal speed of the sound, respectively. 
The wave equation for the deformation parameter Uj (z, t ) is of the form: 


(T Ut _ d 2 o i 
dt 2 dz 2 


(3) 


Taking into account (2), equation |3j for the renormalized deformation, Uj (z, t) looks as follows: 

1 d 2 Uj (z, t) _ d 2 Uj (z, t) 0^ d 2 Nf U, t) d 2 Ej (z) 

c? dt 2 dz 2 pjc 2 dz 2 dz 2 

l r 1 l 


The equation for the defect concentration (interstitial atoms and vacancies) is of the form (7J: 


dN [ i\z,t) 

dt 


— Dj- 


d 2 Nj ] (z, t ) 


dz 2 


-D, 




knT dz 


N"{z,t) 


dUj[z, t ) 
dz 


+ G«- 


A*') ’ 


(S) 


where Dj is the diffusion coefficient of point defects in the z'-th layer, is the generation rate of the 
defects, t^ is the lifetime of the defects in the z'-th layer that is determined by the frequency and the 
amplitude of mechanical fluctuations in the megahertz range (a> 3= 10 6 Hz, r^ ~ 1 ps) that arise in the 
process of the formation of heteroboundaries in stressed nanoheterostructures and in the process of the 
occurrence of defects (acoustic emission) (8). 

As a result, a self-assembled system of equations (4j, i[sj is received for determination of the spatial- 
temporal distribution of the concentration of defects Nv (z, t) and the deformation parameter Uj (z, t ) in 
the different regions of the three-layer nanoheterostructure. 
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The defect concentration can be written in the form: 


N a^» = N do + N >’»’ 


(6) 


where N^J (z, t ) is the spatially inhomogeneous component of the defect concentration. 

Taking into account the presentation jij in the approximation N^J « the equation of diffusion 
is written as follows: 


dN l JJ(z,t) d 2 AT®(z,f) 


dt 




dl' 


dz 2 


D at (0 e d d 2 Ui(z,t ) m 

~ DlN dOk^ dz 2 +G d - 


n 2m. 

Di) 


(7) 


where G^ l) - G® -JvS(z, t)/r® is the generation rate of point defects under the effect of mechani¬ 
cal fluctuations (w ~ 10 6 Hz) that arise in the process of the formation of heteroboundaries in stressed 
nanoheterostructures. 


A further solution of self-assembled systems of equations will be searched in the approxima¬ 


tion 


namely 


2D. 


d 2 Ui (z, t) 
dt 2 


« £o> 


d 2 Uj{z,t) 

dt 2 


= 0 , 


( 8 ) 


d 2 Ui (z, t) 


where L^i is the diffusion length of the defect in the i-th layer. 

In approximation jsj, from equation < 4 J, there will be found 1 ——3 
tion jTJ. As a result, differential equation for determination of the spatial-temporal distribution of defects 
in the stressed heterosystem is received 


and it will be put into the equa- 


dt 


Di 


1 - 


N do 


N 


(i) 


dc 


d 2 N"{z,t) 


dz 2 


- D; 


N dO d 2 £j{z) 
iV®AQ dz 2 

dc 


+ G'® - 


Di) ’ 


(9) 


where iV® — k^Tpc 2 l8^ ] is the critical defect concentration, which being exceeded results in the self¬ 
organization of the defects @. 

In addition, the conditions of the equality of concentration of impurities and their fluxes must be 
satisfied on the boundary layer of the heterostructure shown in figure]]] 


dN%[z,t) 


dz 


= 0 , 


dN^z,t) 


dz 


= 0 , 


MVt-a. t) = N%{-a, f), 


where 


J 2 [a, t) = h {a, t), 

Ji (z, t)--fz |D/iV^fz, f)[l -e,-(z)] |. At the primary moment of time 


'dl 1 ^’ L1 ~ ly dl' 

J\ {-a, t) = J 2 (-a, t), 


( 10 ) 


7V®(z,0) = 0. 


( 11 ) 


Entering the following dimensionless variables 


8 = 


-( 2 )’ 


r (0 _ 




(z, t) = Y t (z, f) G'® t®, 


_ z 
z - —, 
L 


~ L 
L - 

a 


N m = G 'M> 
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P = d 2 


N 


( 2 ) 


dO 


2e 0 


IV® AO V*a 2 Gf ’ 


equations {9} take the form: 

07i(5,0) D 1 


dB D 2 
d Y 2 (5, 6) 


UPY < 

d 
L 
\ 2 


1- 


1V™\0 2 7i(S,0) 


N 


(l) 


dc 


dz 


-Ai 7i(z,0)-l , 


(iPY t 

d 


L 


dB 

073(5,0) L>3 


1- 


N dc) df 1 J 


Il 12) V i 

d 


dB 


D 2 


N^\d 2 Y 3 (z,8) 


1- 


N? 

dc 


dz 


A 3 7 3 (z,0)-1 , 


( 12 ) 


(13) 


where 


|l®/L® j , A 2 = 1 , A 3 = ^ \l} 2) /L® j , and boundary conditions jlO 1 can be written: 


<37i(z,0)i 

dz 


= 0, 


g ( 2 ) t ( 2 ) 

T d 

G (3) t (3) 

(i T rf 


£1 

d 2 


(l [2) \ 2 

d 


A^\07i(-a,0) 


-1 + 


dO 

IV® 

dc 


dz 


d 
L 


N [2 hdY 2 (-a,B ) 


-1 + 


AT® 

dc 


dz 


■ /3a, 


( L( ?) 

2 

( AT® \ 

L 


l Kc) 


5(«,0) D 3 

- 1 - Ba- — 

0z 


V ^ 


^] 07 3 (S,0) 
05 ’ 


at® I 

dc 7 


(14) 


07 3 (z,0)i 
— „ ~ = 0 . 
dz Z—>L 

As seen from equation {12}, parameter /3 describes the nature of the deformation effect caused both 
by the action of the stressed heteroboundary and the action of point defects of the type of compression 
or tension centers. This parameter can take both the positive values /3 > 0 (eo > 0, AO® > 0; £0 < 0, 
AO® < 0) and the negative values ft < 0 (eo > 0, AO® < 0; £0 < 0, AO® > 0). 

The solution of equations (l~3) with boundary conditions ( fl4) is searched in the form: 

7(5,0) = e _Al ' e Z;(S,0), (15) 


where Z,-(z,0) satisfy the following equations: 


dZ\ (z 
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(L 
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2 e 
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i 2 
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dB 


d 2 

1 
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df 
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e Q8+i), 


+ A 3 e A30 , 


(16) 
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with boundary conditions: 


aZi(z,0)i 


dz 


, = 0, 


I z —* — L 

r {2)(2) 

tjj T, 


= -j^e- e ZA-a,e), 

G d T d 
G (3) t (3) 

e- e Z 2 &,0)=-^JLe-WZ3&,e), 

G d T d 


£i 

D 2 


(L^\ 


-Xi6 


-1 + 


K 

N' 


(U\ 

dO 


(1) 
dc I 


dZi{-a,6) 

dz 


( t ( 2 )\ 2 


N 


( 2 ) 


-1 + 


d 0 


N 


( 2 ) 


t r(2) \ 2 


-0 


_ i+ A||a^) 5 = d, 

a* D 2 




-A 3 e 


dc 


-1 + 


dZ2(-a,0) 


dz 


- Pa, 


N d< o 1 5Z 3 (fl, f) 


Nf 

dc 


dz 


(17) 


dZ 3 (z,0)| 


= 0. 


dz 'z —L 

Solutions of equations {16} with boundary conditions {17} are presented in the appendix at 


A i = \ 2 = X 2 = 1, 


r (2)_(2) 
G d T d 

G d T d 


= 1, 


G d^d 

G d V d 


1, * 


2 

d 3 

= 1, — 


D 2 

i a 

D 2 

i a 


= 1. 


3. Analysis of the numerical results and discussion 

In figures [ 2 ] and [i} there is shown the spatial-temporal redistribution of vacancies 7(z,0) (figure [ 2 } 
and interstitial atoms (figure [3} in a three-layer stressed nanoheterosystem GaAs/In, Ga 1 _*As/GaAs un¬ 
der the effect of the deformation caused by both the mismatch between parameters of the contacting 
lattices (£ 0 = A ala = 7%) and by the action of a point defect. Calculations were carried out for the fol¬ 
lowing values of the parameters: eg = 0.07; a = 0.05L, 0.1L (L — the thickness of nanoheterostructure); 


N do ,N dc = °- 5 ; N dl IN dl = °- 8 ; N do IN dc = °- 6 ; c i? = °- 833 q 2 
->(i) _, 


A2) _ , 


( 2 ) _ 


= 0.571 Mbar; T = 300 K; 0 rf 


(0 _ 


5 eV I 


|; Di - 10 5 cm 2 /s; r® = 1 ps. 


0.453 Mbar; = 1.223 Mbar; 


Y[z,e) 



Figure 2. (Color online) Profile of the spatial-temporal distribution of the vacancies concentration in a 
three-layer stressed nanoheterosystem having inhomogeneous-compressed interlayer iV®/AT® = 0.5; 

N do IN dc =°- 8 ; N dl IN dc =°- 6 ^ - 10 - 2 [ formula & 
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Y(z,0) 



Figure 3. (Color online) Profile of the spatial-temporal distribution of the concentration of intersti¬ 
tial atoms in a three-layer stressed nanoheterosystem having inho mog eneous-compressed interlayer 
AT® /AT® = 0.5; IV® / AT® = 0.8; AT® /Af® = 0.6; /3 = —10.2 [formula 


dO dc 


dO dc 


dO dc 


i mog ei 

Bi¬ 


as shown in figures [2] and [3] the profile of the spatial-temporal distribution of the defect concentra¬ 
tion of the type of compression (vacancies, figure[2j or tension (figure[3J centers in a three-layer stressed 
nanoheterosystem is of a nonmonotonous character. If an internal epitaxial layer undergoes an inho¬ 
mogeneous compression deformation (eq < 0) due to the mismatch between the lattice parameter of the 
contacting epitaxial layers, then a decrease (an increase) of the concentration of vacancies (interstitial 
atoms) in the interlayer of the three-layer nanoheterostructures will be observed. 

If the epitaxial layer undergoes the tension deformation due to a mismatch between the lattice pa¬ 
rameter of the epitaxial layer and the substrate ( a s > oq, £o > 0, where a s is the lattice parameter of the 
substrate; ao is the lattice parameter of the stackable layer), the opposite effect will be observed: near the 
heteroboundary there will be accumulation of vacancies and a decrease of the concentration of intersti¬ 
tial atoms. This, in turn, will lead to a decrease of the tension deformation in the epitaxial layer near the 
heteroboundary. 

The effect of impoverishment (enrichment) in the interlayer of vacancies (interstitial atoms) has been 
observed in experimental works (6][10j after the growth (decline) of the intensity of photoluminescence 
in stressed nanoheterostructures. 



Figure 4. The cut of the spatial-temporal distribution of the concentration of vacancies along the growth 
axis at different times: 1 — at the moment t = r®; 2 — t = 5t®; AT®/AT® = 0.5; AT®/AT® = 0.8; 

S d a dO dc dO dc 

]; a = 0.05L. 


23602-7 




















R.M. Peleshchak, N.Ya. Kulyk, M.V. Doroshenko 



Figure 5. The cut of the spatial-temporal distribution of the concentration of interstitial atoms along the 
growth axis at different times: 1 — at the moment t = t®; 2 — t = 5t®; IV^J /A/^ 1 ' = 0.5; = 0.8; 

AT® /AT™ = 0.6; p = 10.2 [formula jlij]; a = 0.1L. 


In figures [3] and [5] numerical calculations of the cut of the spatial-temporal distribution of the con¬ 
centration of vacancies along the growth axis (OX) of heterosystem at different times are presented: 
t = 0 ; t®; 5t^’ (t® is the average time of finding the defect in one of the equilibrium positions in the 
interlayer nanoheterosystem, namely a settled life) and for the different thicknesses of the interlayer of 
a nanoheterostructure (a — 0.05L, figure[4]and a = 0.1L, figure[5j. 

As seen from figures [ 4 ] and |sj during the time interval 0? ts 5 t® there occurs a spatial-temporal 
redistribution of the defects, so that in the inhomogeneous-compressed interlayer [figure [I] region 2, 
formula 111] they become smaller relative to their initial average value Al|S by ~ 13.7%, 16% at different 
times Tj,5Tj\ respectively. Starting from the time t > 5r®, there practically establishes a stationary 
state of the distribution of the defects in a three-layer stressed nanoheterosystem. Thus, the deformation 
field of the interlayer (-a^zsa) clears away the workspace from the defects which finally makes the 
material of the workspace having a greater intensity of photoluminescence flOl . In exterlayers (figure [l] 
regions 1, 3) of the heterosystem, the concentration of the defects asymmetrically monotonously increases 
from the boundary of the contacting materials and becomes larger than their average value AT™, N^j. 

If there are no defects (N^ = 0) in the contacting materials, inhomogeneous deformation is created 
only due to a mismatch between the lattice parameters of contacting materials [e,- (z) = Eq^J = 2], 
and in the absence of the mismatch between the lattice parameters (co = 0), the deformation [f/ w = 
(i 9^ I is caused only by the point defects. 


4. Conclusions 

• It has been established that the concentration profile of point defects at is of non- 

monotonous character with a minimum in the middle of the interlayer In x Gai_ x As which is im¬ 
poverished by the point defect of the type of compression centers when the interlayer of the het¬ 
erostructure GaAs/In x Gai_ x As/GaAs undergoes an inhomogeneous compression, while in the case 
of inhomogeneous tension the opposite effect takes place. Thus, the deformation field and the num¬ 
ber of the defects in the workspace of nanooptoelectronic devices can be controlled. 

• It is shown that if the ratio of the thickness of the middle layer to the thicknesses of the external lay¬ 
ers of a nanoheterosystem GaAs/In x Gai_ x As/GaAs is a/L = 0.05, the ratio of the initial average con¬ 
centrations of the vacancies in the layers of a nanoheterosystem is IN = 0.5, N® /AT® = 0.8, 
N do ^dc ~ 06 ant * va l ue °f the deformation parameter is p = 10.2, then the established con¬ 
centration of vacancies in the middle layer Y(z, 5r®) is less than the initial average concentration 
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by 16 %. If the ratio a/L - 0.1, then the established concentration of the vacancies in the mid¬ 


dle layer is larger than the initial average concentration . Such a reduction of the established 
concentration of the vacancies in the workspace of a nanoheterosystem is correlated with the ex¬ 
perimental results of the work flOl - 


Appendix 

To find the solution of differential equations (16) with boundary conditions (TT) , the integral Laplace 
transformation is used: 


CXJ 

Xi(z, P) = f Zi(z,6)e~ pe d6. 


(18) 


Then, the differential equation (fl6) and boundary conditions (TT) take up the form: 


H\X" [z,p)-pX ] (z, p) + 


P~ 1 


= 0, 


H 2 X” [z, p) - pX 2 (z, p) + = 0, 


P- 1 


H 3 X 2 (z, p) - pX 3 ( z, p) + 


= 0. 


p -1 

Xi(-a,p] = X 2 (-a,p), 

-Hi Xi(-a,p) = -ff 2 X£(-fl,p) ^ 

X 2 {a, p] = X 3 (a, p], 

H 2 X' 2 [a, p) + = -H 3 X' [a, p), 


p -1 


where H, = l-. 


Analytical solutions of differential equations (19} in each layer are as follows: 

1 


(19) 


( 20 ) 


X 1 (z, P ) = C 1 exp|,/-iJ +p(p i) . 

* a p) = C 2 exp I , / JL j) + Cexp ( + X±L. , 
X s fe P ) = C 4e xp + 


-az^-l , 
-a z^a, 

a^z^-l. 


( 21 ) 


Integration constants Q, C 2 , C 3 , C 4 are determined from the following system of algebraic equations: 


Ci exp | - - C^exp y ^ a j ~ c 3 ex P | 


P 


p(p-i) ’ 


-C, /Slexp (- + C 2 /Hxexp (- ^s) - C 3 /Siexp (/jCi) = , 

C2 “ p (\/i s ) +C3exp (- \fk“) - c<exp (- \fk “)=- ^T) ■ 

-C 2 v^exp ( + C 3 v^exp (- ^s) + Q ^exp (- /j|a) = ^71) • 


( 22 ) 
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By carrying out the inverse Laplace transformation we obtain the spatial-temporal redistribution of 
the point defects in the first, second and third stressed layers, respectively: 


Yi(z,0) = e~ e Z\(z,Q) - e“ 0 + 1, 
7 2 (S,0) = Z 2 (5,0)( 1-e" 0 ) (1 + /3), 
Y 3 (z,6) = e~ e Z 3 (z,6) - e -0 + 1, 


(23) 


where functions Zi(z,0), Z 2 (z,0) and Z 3 {z,0 ) are the solutions of differential equations 1 16 with bound¬ 
ary conditions {H}’ respectively: 


where 


Zi(z,0) - - 

b 


e 6 - 1 + e 0 


s/m 


-Erf 


h 


h + £ £ (-D“ 


i= 1 m =0 


fom+1 


X I2i-ir2i-i9 k 1>i(z,6,k)+ X (-D' 


fc=i 


m =0 


fom+1 


X l 2 r 2 9 k ® 2 (z,9,k), 


k= l 


<Pi(z,0,fc) = 




46 




r 2 i-iT (k + 1 ) 


- (- 1 )' 


s/0r 2 i-iT[k+ l) 


4TI l r l 
+ 




0 r(fc+i) 


(24) 


(25) 


<f> 2 (z,0,fc) = -- 




•2.’ 40 


VQT[k+\) 


r 2 r(fc+l) 


a = (y/Hl- / h 2 )( v a h 2 - Vm), b=[Sm+ y/m)(y/ih+ Vm), 
h=ps/m[Vm- \fm}> i 2 = --^=, 

h=pVm(Vm+ y/m). h =^fr 2 ’ 


z + a 4a(m+l) 

ll ~~Vm + mm 


z + a 2a(2m+\) 

r2 ~~^m + s/m 


z+a 4am 

~Vm + 7m' 


(26) 


where 


Z 2 (z,0) = ^ [(e 0 -l)(p 3 + p 4 ) + (p 5 + p 6 )e 0 Erf(v / 0)] 

4 oo sjtn oo 

+ X X 


um+1 

i=lm=0 D jfc=l 


X PieiB k <& 3 {z,e,k) } 


$> 3 {z,d,k) = -- 




’2’ 46 




e,T(fc+ 1) 


1,-^1 iFill-fc,!,-^ 


’2’ 46 


v/#2 v / 0e,T(fc+ g) 


er(it+|) 


(27) 


(28) 
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Spatial-temporal redistribution of point defects ... 


Pi^p7m[mm- 7m), P 2 = pa(^ir 2 - mm), P3- 


P2a 

\fih' 


Pi- 


~ P \fl~h ( \Zl~h + 7m), 


P 5- 


p^a 

7m’ 


P 6 


■ P7m(7m- 7m), 


z + a 4a(m+l) 

ei_ ~7^ + 7m ' 

z-a 4a(m+l) 

^3 = ,_ + -;=-, 

\fih \fih 


z+a 2a{2m+l) 
62 =-== + -;=-, 

mm mm 

'z-a 2a(2m+l) 

64 = —= +- -= — . 

mm mm 


where 


Z 3 (z,0) = 


1 

b 

OO 


e e - 1 + e 0 —^=Erf( 7d) 


7m 


2 00 




i =1 m-0 


fom+1 


L fe*-ir2*-i0*<E > 4(S,0»fc) + £ L kr20 k O 5 {z f 0 t k) 


it=i 


^ m+1 a 


i^Fil iFif-fci,-# 


® 4 U,0,fc) = -- 


Ver[k+l) 


r 2 i-iT{k+l) 


iFi^-k, 1 ,,-^) iFifl-U,-%r 


n /#2 v / 0r 2 /_ir(fc+1) 


0 r(fc+|) 


(29) 


(30) 


< t> 5 {z,e,k) = -- 




’2’ 40 


V0r(fc+i) 


r 2 r(fc+ 1) 


/i = /3v^(v^i- \/*fe), fe = - /lfl 
I 3 — P 7 ~m> ( 7 ~m. + 7 ~m) > (4 = 


7m’ 

I 3 C 1 

7m’ 

z-a 4a(m+l) z-a 2a(2m+l) z-a 4am 

ft = —== H-;=— , r 2 = —— = H-— = —, r 2 = ——= h —= . 

^ ^ mm mm 


mm mm ’ 


mm mm ’ 


(31) 
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npocTopoBO-HacoBi/iM nepepo 3 nofli.n tomkobmx fle<|>eKTiB y 
Tpn lli apoBi/ix Hanpy>i<eHMX HaHoreTepoci/icreMax y Me>i<ax 
caivioy 3 rofl>i<eHoi flec|>opMa^MHO-Ai/jc|>y 3 mHoi MOfle/ii 

P.M. rie/iemaK, H.fl. l<y/n/iK, M.B. flopoLueHKO 

flporo6m4bKMti AepxaBHi/iw neflaroriMHUM yHiBepcmeT iM. iBaHa cppaHKa, 

By a. I. <hpaHKa, 24, 82100/(poro6u4, YitpaiHa 

no6yflOBaHO Moge/ib npocropoBO-HacoBoro po3nofli/iyTOHKOBHxfle(|)eKTiByTpnujapoBiti HanpyxeHiii HaHore- 
TepocucTeMi GaAs/In x Gai_ A -As/GaAs 3 BpaxyBaHHHM caMoy3rofl>KeHoi flet|)opMaL(iwHO-flH(|)y3iCiHoi B3aeMOflii. 
y Mexax u,iei' MOfle/ii po3paxoBaHO npocf)i/ib npocropoBO-HacoBoro po3nofli/iy BaxaHCii) (Mi>KBy3/ioBnx aTOMiB) 

y Hanpyxemi) HaHoreTepoa'icreMi GaAs/IrixGai-^As/GaAs Ta noKa3aHO, mo y BwnaflKy BCTaHOB/iem-ifl CTapio- 

( 2 ) 

HapHoro CTaHy (f > 5 t^ ), KOHLteHTpapiti BaxaHcin b HeoflHopiflHO CTHCHyTOMy BHyipiuiHbOMy Luapi e MeHLLiOK) 

( 2 ) 

BiflHOCHO BHxiflHoro cepe^Hboro 3Ha4eHHti N\‘ Ha 16%. 

KsiiOMOBi c/ioBa: npocTopoBO-nacoBwi/i nepepo3nofli/i, BaKaHCi'f, Mi>KBy3JioBi aroMia 
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